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Scalings of the synchrotron cut-off and turbulent 
correlation of active galactic nucleus jets 
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ABSTRACT 

We propose a new analytic scaling of the cut-off frequency of synchrotron radiation 
from active galactic nucleus (AGN) jets that are nonuniformly filled with many fil- 
aments. The theoretical upper limit is provided independent of magnetic intensity, 
spectral index, coherence and correlation length of filamentary turbulence, etc., such 
that v c ~ 6 x 10 20 <5[(r— l)/r] 4/,3 (&/10 -4 ) Hz, where 8, r and b are the Doppler beaming 
factor, shock-compression ratio and energy-density ratio of the perturbed/local mean 
magnetic field of the filaments, respectively. Combining our results with observational 
data for 18 extragalactic sources, a constraint on the filament correlation length is 
found, in order to give the number scaling of filaments. The results suggest that, 
in particular, the jets of compact BL Lacs possess a large number of filaments with 
transverse size scale smaller than the emission-region size. The novel concept of the 
quantization of flowing plasma is suggested. 

Key words: acceleration of particles - radiation mechanisms: non-thermal - turbu- 
lence - methods: analytical - galaxies: evolution - galaxies: jets. 



1 INTRODUCTION 

The pronounced extension of the synchrotron radio contin- 
uum to the X-ray region, not formerly predicted, is now 
kno wn to be a very common f eature among radio jets 
(e.g. lHarris fc Krawczvns ki 2006). The crucial point that 
should be emphasized is that, in many sources, still no 
clear signs of sharp cut-off are found in the Chandra regime 
(e.g. 3C 66B: [Hardcastl e et al, I l200ll : M 87: iMarshall et al.l 
120021 : 3 C31:lHardcastle et alj|2002l; M 84: lHarris et al.ll2002l : 
3C346: IWorrall fc B irkinshawl lioOot ). and there is even an 
implication of the possible appearance of a gamma-ray tail 
in the synchrotron spectral component. This fact is of con- 
siderable interest, directly related to the challenging issue 
of how high an energy range electrons could be accelerated 
up to in situ. In a conventional fashion, the electron energy 
distribution is truncated at an assumed 'highest energy', 
which could be responsible for the observed spectrum. Inde- 
pendent of such a makeshift solution, however, it is desirous 
to elaborate on a particle acceleration theory that can nat- 
urally account for the extended continuum, with reference 
to the morphological details innovated by the recent very- 
long- baseline interferometry survey (e.g. lLobanov fc Zensua 
l200ll ). From another point of view, theoretical expansion 
would also be requisite for corroborating a unified scheme of 
radio -loud active galactic nuclei (AGNs; lUrrv fc Padovanil 
1 19951 ) . None the less, the related thorough study required 
to elucidate a universal mechanism underlying the many 



appearances of the energetic continuum has not been per- 
formed as yet. 

In a simplistic model assumin g a homogeneous magnetic 
field over a particular system, iBiermann fc Strittmatten 
i|l987l ) first considered a combination of small pitch-angle 
scattering diffusion and Fermi-type acceleration of electrons 
to evaluate the frequency of synchrotron emission from the 
highest energy electron: 



= 3 x 10 14 [3b{U/cf] Hz, 



G) 



where U, c and b (< 1) are the shock speed, the speed 
of light and the ratio of the perturbed to the mean en- 
ergy density of the global magnetic field, respectively. Al- 
though a symptom of cut-off in the range of equatio n |T]| 
has been suggested in some earlier publications (e.g. [ Keell 
ll988l : |Perez -Fournon et al. l ll988l : lMeisenheimer et al.lll996l ). 
the broad-band spectral fitting is now found to, if anything, 
entail a break followed by the aforementioned high-energy 
extension, unlike a sharp cutoff. A clue aiding in achieve- 
ment of the required higher acceleration efficiency might be 
the fine structure like filaments inside jets, which can be rec- 
onciled with th e observational d e tailed morphology. Based 
on this notion, iHonda fc Hondal (|2005al ) have investigated 
the root-mean-square (rms) diffusion of electrons in inho- 
mogeneous magnetic fields, and found that the intrinsic fre- 
quency of the diffuse synchrotron component could reach as 
high as 
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v c ~ 10 24 B(mG)[d(pc)] 2/3 (U/c) i/:i Hz 



(2) 



where B and d are the rms magnetic field strength and trans- 
verse correlation length of filaments, respectively. The rele- 
vant argument focusing o n the specified case of the nearby 
M87 jet has been given (jHonda fc Hon da 2007), although 
the value of d remained unsolved. Its determination is nec- 
essary not only to estimate equation {2} but also complete 
the filamentary jet model, w hich could involve various radi- 
ation channels l|Hondall2008h . 

This paper has been prepared to spell out explicitly 
the scalings of the cut-off and break for the dominant syn- 
chrotron spectrum established via normal (non-diffuse) pro- 
cesses in filamentary jets, and to revise equation (fT)l. includ- 
ing relativistic beaming effects. Making use of the results, 
we attempt to extract the scaling of the lower limit of the 
filament correlation length, d m i n , from observational data 
for sample extragalactic jets. We find the property that the 
rfmin value increases as the size of emission regions increases. 
A corollary derived from this is that the allowable maximum 
number (capacity) of filaments with the outer scale (~ d) in- 
creases as the propagation distance of the jet increases. For 
the situation in which AGN jets carry huge currents driven 
by th e central engines of black holes l|Appl fe Camenzindl 
1992), we take the current filamentation into consideration 
in order to calculate the number of filaments (N). From 
comparing N with the capacity of the outer-scale filaments, 
we infer the population of fine filaments inside the jets. As a 
result, it is demonstrated that compact BLLac objects such 
as Mrk421 and 501 would possess typically ~ 10 11 filaments 
with various transverse size scales that are smaller than the 
entire emission-region size. 



2 THEORETICAL MODEL AND ANALYSIS 
2.1 The filamentary jet model 

According to the original idea in iHondafc Hondal (j2004r i. we 
rely on the hypothesis that a jet consists of many magne- 
tized filamejrts^ccxHnmodated by some radio observations 
(e.g. lOwen et alJll989h . For information, the circumstance 
concerned is illustrated in Fig. [1] The quasi-stationary mag- 
netic fields are considered to be generated by leptonic cur- 
rents, and retain an energy-density level co mparable to t hat 
of the random component of lepton motion l|Hondall2007l ). It 
is then expected that the Poynting flux level is less than (or, 
at most, comparable to) the lepton energy flux. In respect of 
the global energy budget, the small portion that takes part 
in controlling the transverse dynamics (Section [5]) is loaded 
from the energy reservoir of the jet bulk. If the dominant 
kinetic energy is carried by hadronic components then an 
excess of lepton energy departing from the mass ratio (say, 
0.05 per cent for a simple electron-proton plasma with no 
pairs) is entailed in order to generate the currents and mag- 
netic fields, in acc ord with arguments conce ring the power 
of blazar jets (e.g. ICelotti fc Ghisellmill200sh . 

The merging of current filaments can be associated 
with the inverse cascade of turbulent magnetic energy. For 
the inertial range of the turbulence, we invoke the phe- 
nomenological expression for local magnetic intensity B = 
-B m (A/d)" 3_1 ^ ,/2 , where A (^ d) is the transverse size of a 
filament (Fig. [1} and /3 corresponds to the turbulent spec- 
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Figure 1. Schematic of the AGN jet that contains bright knot- 
like features associated with shocks. The magnification stands 
for a local sample around the working surface accompanied by 
electron acceleration and synchrotron emission. 



tral index. As seems plausible, this (zeroth order) mag- 
netic field, trapping lower-energy electrons, is disturbed by 
Alfvenic waves to allow resonant scattering diffusion. Also, 
in analogy to the knot and hotspot features, it is reasonable 
to suppose a relativistic (non-relativistic) shock overtaking 
the relativistic (non-relativistic) flow with the Lorentz factor 
Tj = (1 — /3 2 ) -1 '' 2 , such that the shock viewed in the up- 
stream rest frame likely has a non-/weakly relativistic speed. 
Providing the fields defined upstream, we consider the stan- 
dard diffusive accel eration of bound elect rons due to the 
Fermi-I mechanism IjHonda fc Hondall2007h . which yields a 
power-law energy distribution (oc 7 _p for 7 ^ 7*). 



2.2 Updated scalings of break and cut-off 
frequencies 

The acceleration time-scale t aC c is of the order of the cycle 
time for one back-and-forth d ivided by the energy gain per 
encounter with the shock (e.g. Gaisscr 1990). We here adopt 
a standard expression for t acc that involves a negligible con- 
tribution f rom particle escape downstream to the cy cle time 
and gain (|Ostrowskilll988l ; iHonda fc Hondall2005bh . What 
we deal with below is, therefore, the non-ageing (or flaring) 
regime around the working surface (a related issue is revis- 
ited in Section 13. ip . Compatible with this, it is reasonable 
to suppose that, for major electrons bound to the outer- 
scale filament with maximum field strength ~ B m , the ac- 
celeration is degra ded by strong synchrotron cooling rather 
than escape losses i|Hondall200cf ). The synchrotron emission, 
largely polarized, appears to constitute the dominant radio 
continuum with a power index of a = (p — l)/2 in the flux- 
density spectru m (S v oc v~ a , provided the electron density 
is uniform, e.g. lLongairl fl994). The spectral break reflects 
the highest energy of an accelerated electron ( r y*\\^dmc 2 ) 
at which f acc is comparable to radiative time-scale. From 
balancing these time-scales, we can analytically derive the 
maximum Lorentz factor, 76 = 7* \\ = d (along with the guide- 
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lines provided in lHondall2008h ; for convenience, the generic 
expression is explicitly written as follows: 

lh = [f(b,r;B m ,d)]^ 3 -^ [g(B m ,dT^'^-^ . (3) 

Here, b is the energy density ratio of the perturbed/local 
mean magnetic field of filaments, which is assumed to be a 
constant smaller than unity: 1 (checked below), r is the 
shock-compression ratio, /3'(> f) is the turbulent spectral 
index of the Alfvenic fluctuations (superimposed on the local 
mean field) and the dimensionless functions are defined as 



f = 9n 2 (/3'-l)b- 



1 B m d 2 



eB m d 



(4) 



e ' 2mc 2 ' 

where c, e and m are the speed of light, elementary 
charge and electron rest mass, respectively. The ob- 
served break frequency can be simply estimated as vt = 
(3/47r)^[eB TO /(mc)], where 5 = rj\l - ft cosfl) -1 is the 
beaming factor. When setting the standard value of ft' = | 
(for Kolmogorov turbulence), we obtain the scalings 



7t , = 1.1 x lO 6 ^ 4 ^ 



5/4.-1/2 

Li , 



(5) 
(6) 



v h = 5.2 x 10 15 S^B^^d" 1 Hz, 

where £_ 5 = £/l(T 5 and £ = 6[4 (r - 1) /3r], B m in mG and 
d in pc. Note that both equations |5| and © are indepen- 
dent of p. 

In filaments smaller than d, the bound electrons can be 
accelerated up to a higher energy, i.e. 7*|A<d > 7b, owing 
to the weaker synchrotron loss. While the synchrotron flux 
density is lower, the power-law tail is retained up to higher 
frequency (reflecting the enhanced 7* value). Apparently, 
this property can be responsible for the extended contin- 
uum going beyond u^. However, the increase in acceleration 
efficiency ought to be limited at a critical coherent length 
A c , below which the transverse escape of electrons dominates 
the radiative loss. In the escape-dominant regime, another 
7* scaling is derived from the spatial limit condition, i.e. 
equating the electron gyroradius r g with a(A/2), where a is 
a dimensionless factor smaller than unity (see Section 12.31) . 
The equation for 7*|a<a c = 7*|a>a c contains the solution of 
A = A c at which 7* takes a peak value, to give the expression 
Ac/d = a 2(3-£')/(3/3+U / a/(s/ > +i) ff -8/(s/ > +i) < Substituting the 

Ac expression into the 7* scaling, the critical Lorentz factor, 
defined as 7 C = 7*|a=a c , is obtained as 



7c = a 1 



[(/3+l)/3'-2]/(3/3 + l) j(/3 + l)/(3/3+l) 5 -(/3+3)/(3/3 + l) 



(7) 

and, in turn, the intrinsic cut-off frequency (in 
the observer frame) can be estimated as v c = 
(3/47r)57 2 [eB m /(mc)](Ac/d) (,9 " 1)/2 , which is found to 
be recast in the simple form 



v c = 277r(ft - l)Sa fi 



1 mc* 



(8) 



independent of ft B m , d, etc. For ft = 5/3, the quantities 
scale as 

7 c = 3.8 x 10 s a^e^B^d 1 ^ 7 (9) 
for = 2 (see lMontgomerv fc Liulll979l as an example) and 



= 4.5 x 10 19 a 2/3 5£_ 5 Hz 



(10) 



respectively. Note that equation (|10|) is independent of ft In 
the ordinary case in which the diffuse emission level is lower 



(cf. the discussion in lHonda fc Hondal |2007i), equation (|10[) 
(instead of equation [2j provides the observed cut-off fre- 
quency of the synchrotron radiation, which is of non-diffuse 
at least in the scale of ~ A c . Note that the absorption in 
extragalactic background lights is expected to be insignifi- 
cant, as long as the radiation freque ncy is in the range below 
~ 10 24 Hz (e.g. iKneiske et alj|2004h . 



2.3 The smearing effect on the spectrum at v c 



The large v c value expected in equation (|10[) is owed to the 
efficient acceleration of electrons bound to the local mag- 
netic fields of filaments with size-scale ~ A c . In this regime, 
the time-scale of escape from a small-scale filament com- 
petes with the radiative loss time-scale, and therefore we 
provide additional discussions concerning the spatiotempo- 
ral property of diffusion loss, in order to identify the value 
of a. The escape is dominated by diffusion across the mag- 
netic field line permeating through the current filaments, 
so that for the electrons confined within radial size ~ A/2, 
the escape time is estimated as t CBC ~ (X/2) 2 /c£±, where 
£± is the mean free path (mfp) for diffusion perpendicu- 
lar to the magnetic field line (e.g. Hillas 1984). Recall that 
the perpendicular diffusion is characterized by the coeffi- 
cient k± — K||/(l + rj 2 ), where r\ — l\\/r g , fty = \r)r g c 
and t\\ = r g /[b(P' - l)][A/(2r 9 )]' 8 '- 1 (for > 1) denotes 
the mfp for the diffusion parallel to the field l ine (e.g. 
iMiicke fc Protherod l200ll : iHonda fc Hondal l2005bh . Along 
with this, we invoke the relation £± ~ V~ 2 ^\\ f° r V 2 ^ 1> 
which is the region of interest to us, and then the escape 
time can be expressed as t CBC — (£^/c)(\/2r g ) 2 . As a result, 
we find that equating t eS c to t aC c correctly provides the afore- 
mentioned spatial limit equation in the form of r g = a(A/2), 
where 



y/{r-l)/r. 



(11) 



One can ascertain that equation (|lip is also valid for /?' = 1 
(not shown). In the strong-shock limit we have a = 0.87, 
and suggest that setting a to unity is nearly adequate for 
the simple treatment of escape loss. In the case in which 
the mfp related to the diffusion considered is anomalously 
longer than the estimated range of £±, the a value becomes 
effectively smaller than that of equation (|11[) : i.e. 7 C and u c 
could reduce, even if the inverse Compton (IC) and redshift 
effects considered later are negligible. The uncertainty in a 
is expected to cooperate to smear out the actual spectrum 
around u c . 



2.4 The inverse Comptonization effect 

There exists another possible effect that prevents electrons 
from being energized up to the range given in equation Q. 
The starlight and dust emission that emanate from host 
galaxies potentially serve as external targets for the in- 
verse Comptonization of accelerated electrons. For rela- 
tivistic jets, one may consider these and also the cosmic 
microwave background, because the radiation energy den- 
sity is boosted in t he comoving frame by a factor of V 2 
ijStawarz et all [2003). Particularly in the weaker magnetic 
intensity regions of A < d, the comoving radiation energy 
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density (denoted as ii ra d) could be comparable to the mag- 
netic energy density, so that the time-scale for Thomson 
scattering loss could compete with that for synchrotron loss. 
It is pointed out that similar circumstances can appear even 
for a synchrotron originator, if the non-local effect is sig- 
nificant for radiative transfer. On the other hand, it is un- 
likely that inverse Comptonization affects the lepton ener- 
gization to the 7 C level of equation Q, because the Klein- 
Nishi na effect tends to redu ce the scattering cross-section 
(e.g., iHonda fc Hondal I2007T I . In any case, the competition 
between these conceivable radiative channels is complicated, 
dependent on the detailed parameters inherent to sources. 

With this aspect in mind, we here simply clarify the 
maximum level of IC losses that can affect the electron ac- 
celeration up to 7 C . For example, one can identify the max- 
imum equivalent IC field (denoted as B lc ) corresponding to 
the quantity \/87ru ra d by conservatively taking the balance 
of t acc with the Thomson time-scale at j c . Considering the 
synchrotron limit at jb, it is found that any IC losses can 
safely be ignored (i.e. equations [9l and 1101 are valid) if 



fin 



< 



(/3-l)(3-/3')/[(/3 + l)/3'-2] 



(12) 



The criterion for {fi,P') = (2,5/3) reduces to B ic /B m < 
(ib/lc) 4 ^ 9 , which would generally be satisfied. In the follow- 
ing, therefore, we are concerned with the negligible contri- 
bution of Compton losses to the maximum energy analysis, 
aiming at providing a generic theoretical upper limit on the 
(non-diffuse) synchrotron radiation frequency that has not 
been manifested so far. 



3 COMPARISON WITH OBSERVATIONAL 
DATA 

3.1 The ageing effect and degenerative v c redshift 

We further expand our discussions on the validity of the 
scalings of Vb and u c , in light of the comparison with obser- 
vational data. In equation (|6]) , we have the Vb scaling appar- 
ently proportional to B m 3 ^ 2 . However, there is a potential 
that the correlation length-scale d is virtually related to B m , 
yielding an explicit magnetic intensity dependence of Vb- To 
see this, in Fig. [2] we plot magnetic field strength B against 
emission-region size D for 48 sample objects (102 features) 
that include blazars, Fanaroff-Riley type I /II radio galaxies 
(denoted as FRI/II, respectively) and quasi-stellar objects 
(QSOs). It is found that the power-law fitting suggests a 
generic scaling of D ~ £j-2.i~(-2)^ though the translation 
to a d — B m relation is not trivial. At this juncture, provided 
the scaling is the same, one may read off Vb ~ fim 5 ~°' 6 j 
when supposing a simple ordering of B/B m < 1 and that 
d ~ D ' 7 is retained (cf. Section [4|, we preliminarily ob- 
tain Vb ~ S~ 01 ~ +0 ' 1 . In any case, the resulting magnetic 
intensity dependence of Vb is thought of as being weak in 
the present context. This property appears to be, if it is 
corre c t, markedly differen t from Vb oc B~ 3 (|Brunetti et al.l 
120031 ; ICheung et all 120051 ), which was derived considering 
the synchrotron cooling of ageing electrons in hot spots 
jMeisenheimer et al.lll997l ). We point out that observations 
provide some support for the classical scaling (for hotspots), 
implying that the present analytic results cannot directly be 
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Figure 2. Magnetic field strength B, versus emission-region size 
D, for blazars (filled diamonds), FRI (filled squares) and FRII 
(open squares) radio sources and QSOs (crosses). B equipar- 
tition values in the rest frame, i.e. B cq s—i8~ 5 / 7 , are from 
IKataoka fc Stawarj l |20 05j) (KS05), except for the marginal data 
of four compact BLLacs, which are commonly from the syn- 
chrotron self-Compton models provided in the literatures in Ta- 
ble [T] It is noted that for the 'SYN' features classified in KS05, B 
is simply set to B eq s=l- The solid line indicates the best-fitting 
scaling B = 1.7D-°' 48 . 



applied to hotspot phenomena, although the applicability of 
the filament model as such still remains unclear. 

In the synchrotron spectrum, the appearance of ageing 
effects might be pronounced, particularly during the low ac- 
tivity around the working surface, which can be translated as 
the quenching of turbulent energy injection at the smallest 
size-scale (i.e. the largest wavenumber, fc m , i n the filamen- 
tary turbulent spectrum: IHonda et all feoOO). The inertial 
inverse cascade will lower the upper cut-off k m , resulting in a 
lack of smaller-scale filaments and magnetic energy conden- 
sation around the outer scale ~ d. When entering into the 
degenerate regime of k^ 1 > \ c /2ir, equation (|10|l is violated, 
and the further cascade leads to the 'redshift' of cut-off fre- 
quency moving toward the lower limit compared with the Vb 
range given in equation ((6]l . This scenario provides a reason- 
able interpretation of the observational fact that there are 
still many obj ects with a cut-off signature around the opt i- 
cal band (e.g. iMeisenheimer et~ai]|l997l ; iMack et all 120091 ") . 
while there is evidence of sy nchrotron X-ray emission in a 
number of jets and hotspots (Ha rris fc Krawczvnskil feOQG); 
i.e. the inspected low-energy cut-offs do not always imply 
extremely small values of £ in equation (|10p . In the specific 
regime, a spectral break below the low-energy band could 
reflect the oldest electron population along the classical sce- 
nario, whereupon the scaling of Vb oc B~ 3 might be realized 
l|Brunetti et al.ll2003l ; ICheung et al]|2005l ). 



3.2 Theoretical upper limit of v c 

Listed in Table Q] are the relevant data for radio sources, 
which are used to evaluate the value of v c in equation (|10|l 
and to compare the measured v\, with equation © . Basically, 
we have sought X-ray sources that retain the extension of 
the synchrotron continuum above the confirmed Vb, although 
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for the data selection the statistical bias might inevitably 
be involved. The samples selected here include BL Lac ob- 
jects, FRI radio galaxies, core-dominated quasars (CDQ) 
and a gamma-ray quasar (GRQ); these are found, in par- 
ticular, to provide a well-defined parameter set that could 
also be referred to the correlation length estimate later in 
equation (fl3|) . Taking a conservative approach, FR1I sources 
have been excluded at the moment, since hotspot phenom- 
ena in the tips often incur the ageing effects that have been 
annotated in Section [3.f I and also sho w some peculiar fea - 
tures in energetic emission (e.g. PicA; lTingav et al.l l2008). 
On the other hand, there are som e cases for which FR II fea- 
tures resemble knot s in jet (e.g. iMeisenheimer et all 1 19971 ; 
iMarshall etai]|2010h . and therefore the applicability of the 
present model to FRII sources will be investigated in more 
detail, though it seems to be somewhat beyond the scope of 
this paper. 

The updated X-ray data provide the lower limits on v c 
in equation (|10[) . and thereby the lower limits of the unfixed 
parameter, £ m i n , for given a ~ 1 (upper limit) and measured 
8 value s. Then, by invokin g a simple relation p ~ (r + 2)/(r — 
1) fsee ISchlickeiserl 12002. for details), the lower limits of ft 
can be estimated such that 6 m i n ~ [(p + 2)/4]£ m j n . 

Fig. [3] plots the predicted values of v c against ft (in the 
range of ftmin), except for CDQs (for reasoning, see Sec- 
tion |4| . According to equation |TT}, the value of a is sim- 
ply set to y/3/(p + 2) hereafter. It is confirmed that the 
condition of b < 1 (introduced in Section 12. 2|) is satis- 



fied as long as v c sC 10 Hz. For a given b, the v c val- 
ues of the compact BL Lacs are likely higher than those 
of the FRI sources, mainly because of the higher 5 val- 
ues for the B L Lacs. For the set (p. 6) = (1.6, 10) (Mrk 421 



Mrk501: e.g. iTayecchio et al.l l200ll : iKonopelko et alj|2003l; 



iBIazeiowski et aT 120051) . a specified value of ft = 1 



x 10" 



yields u c = 4.7 x 10 Hz, in conformity with the extremely 
high-frequency ranges of synchrotron cut-off expected in 
TeV gamma-ray emitters (|Hondall2008h . It is also mentioned 
that the scaling of a beaming GRQ (P KS 0208-512 with a 
hundred- kpc scale; ISchwartz et al.ll2006f) appears to be c lose 
to that of a kpc-scale BL Lac (3C 371: |pesce et al.ll200ll . not 
shown in the figure). 



4 CORRELATION LENGTH-SCALE AND 
'PACKAGE' OF FILAMENTS 

A worthwhile manipulation is to combine equation © with 
pup , eliminating the parameter £. We solve this for the cor- 
relation length d to find the scaling 



d = 5.4 x 10" 



-1 r-l/2 

d ' 1 



3/2 d-3/2 

c ,17 m 



pc, 



(13) 



where 1^14 = ^(,/10 14 Hz and v c> vr — i^ c /10 17 Hz. By putting 
the parameter values of Table [T] into the right-hand side of 
equation (|13[1 . one can fix the allowed d domain for each 
sources. In the derivation of equation (|13|l . it is supposed 
that the param eter £ in equa tions ([6]) and (|10l) takes a com- 
mon value (cf. iHondal [20081 ) . If the back-reacti on effect of 
accelerated particles on the shock structure fe.g. lBlaslll2002l ; 
iKang fc Jones) 120071) comes into play with a positional de- 
pendence, the value of £ may be different in between the 
regions of A ~ d and A c , as ft and r change spatially. Related 
to this issue, the strong non-linear effect also invalidates the 
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Figure 3. The cut-off frequency v c of synchrotron emission from 
the extragalactic jets (labelled) as a function of the energy den- 
sity ratio of the perturbed/local mean magnetic field of filaments, 
b. The darkly and lightly shaded areas indicate the expected do- 
mains for Mrk 421 and 501 with (p, 6) = (1.6 - 1.8, 10 - 100) and 
M87 with (2.3, 1 - 3), respectively. 



assumption of 6 < 1 followed by r\ 1. We mention here 
that the back-reaction effect may scatter the results shown 
below, resulting in an effective lower estimate of d (in equa- 
tion I13[) , but quasi-linearity (6 <C 1) seems to be almost 
satisfied even in the smaller filament with A c , at which the 
acceleration efficiency is maximum. 

Fig. |3£a) plots the lower limit, d m i n , against the size 
of emission region D. Note the relation d m in < D, which 
maintains the theoretical consistency. As seen in the figure, 
the plotted points including bars may be separated roughly 
into three groups (labelled as BLLac, FRI and CDQ). No- 
tice that the FRI group apparently absorbs a large-scale 
BL Lac (3C 371), implying the similarity of turbulent states. 
Although the uncertainty ascribed to the observations is not 
small for the moment, it seems natural to suppose that the 
(D,d m in) points of the BL Lac and FRI share a common 
power-law scaling. The CDQ group is apparently isolated in 
a marginal region below the scaling (mainly because of the 
larger uncertainty of v c , stemming from the flat spectral fea- 
tures; cf. Table [TJ, although the d m i n levels are still retained 
above a skin depth for the plasma density of > 10~ 6 cm -3 , 
satisfying the restriction A/d m j n ^ 1 (cons istent with the 
filamentation instability; iHonda et aDl2000h . With the ex- 
ception of the marginal data, the power-law fitting yields a 
scaling of d m i n , P c ~ 5 x 10 _5 Dp C 7 , where the units are pc. 
This arguably reflects the history of how the filamentary tur- 
bulence evolves as the jet propagates, increasing the width. 
The physical implication can be revealed by translating a 
quantity < D 2 /d^ nin as the upper limit of the number of 
outer-scale filaments, (N\^d)max- 
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Table 1. Sample parameters to evaluate equations l|10[ l and 1 3 fl - 



Host/feature (class) 


D(pc£ 


Pi 


Bm(mG£ 


S 


v b (Hz) 


y c (Hz) 


References^ 


1ES 1959+650 (BL Lac) 


3.8 x 10~ 3 


2 


400 


20 


10 16 - 10 17 


> 10 19 


1 


Mrk421 (BL Lac) 


3.1 x 10~ 3 


1.75 


1000 


55 


10 16 


> 10 19 


2,3 


Mrk501 (BL Lac) 


2.3 x 10~ 3 


1.60 


400 


50 


10 16 


> 10 20 


2,1 


PKS 2155-304 (BL Lac) 


1.5 x 10~ 2 


1.35 


1400 


28 


(.5 - 1.3) x 10 16 


> 10 la 


5,6,7,8 


3C 371/A (BL Lac) 


252 


2.52 


.81 


6 


10 14 - 10 lj 


> 10 17 


9,10,11 


M84/N3.3 (FRI) 


8 


2.3 


1.3 


.6 — 1.25 


4 x 10 13 


> 2 x 10 17 


10,11,12 


3C31 (FRI) 


68 


2.1 


1.1 


1.3 


10 13 


>2x 10 17 


10,11,13 


3C 66B/B (FRI) 


85 


2.2 


1.2 


1.2 - 2.8 


10 14 


> 2 x 10 17 


11,14,15 


M87/A (FRI) 


110 


2.34 


5.1 


1-3 


(.27 - 7.0) x 10 16 


> 2 x 10 17 


16,17,18,19,20 


Cen A (FRI) 


154 


2 


6.4 


4 


10 12 - 10 14 


> 2 x 10 17 


10,11,21,22 


3C 346 (FRI) 


320 


2 


7.2 


< 3 


10 12 - 10 13 


>2x 10 17 


11,23 


PKS 0208-512/R2 (GRQ) 


5240 


2.4 


.135 


7.5 


s£ 10 14 


> 10 17 


11,24,25 


3C 273/D (CDQ) 


3000 


2.6 


220 


4.2 


10 13 - 10 14 


> 10 15 


26,27 


4C 49.22/B (CDQ) 


7800 


2.4 


.75 


14 


10 13 - 10 14 


> 10 15 


10,11,28 


4C19.44/A (CDQ) 


10000 


2.4 


.86 


14 


10 13 


> 10 15 


10,11,28 


PKS 1202-262/R2 (CDQ) 


11400 


2.4 


120 


12 


10 11 - 10 14 


> 10 15 


11,25 


PKS 1030-357/R2 (CDQ) 


13600 


2.4 


223 


9.2 


10" - io 14 


> 10 15 


11,25 


PKS 0637-752 (CDQ) 


16700 


2.62 


2 


10 


< 3 x 10 12 


^ 5 x 10 14 


11,29,30 



a Values set simply to D = L/50, if not well defined. 

" Values inferred from comparing the measured radio spectral indices with a = (p — l)/2. 

c For an ad hoc manner, set to ten times the field strength suggested in Refs (for the reasoning, cf. lHonda|[200sT) . 
d Re f: 1 iKrawczvnski et al. |l 2004l. 2. iKonopelko et"al]l2003l. 3. iBIazeiowski et al.ll2005l. 4 iTavecchio et al.ll200ll. 5. lUrrv et al l 
Il997l. 6. iKataoka et aljbqool. 7. IZhang et al.ll2002l. 8. iDominici et alj 2004 9. iPesce et al-lboOll. 10 IKataoka k. Stawarzll2005l. 



ll.lHarris fc Krawczvnskill2006l. 12.lHarris et al.ll2002l . 13. lHardcastle et al.ll2002l, 14.lTanslev et al.ll2000l 15. lHardcastle et all 



l200ll. 16.IOwen et al.lll989l , 17. IBiretta et al lll99ll. 18. |Perlman et al.ll200ll. 19. IMarshall et alj|2002l, 20lWaters fc ZepJ^OOa 



21. IKataoka et al.ll2006l. 22 lHardcastle et al.ll2006l. 23.IWorrall fc Bir kinshaw 2005, 24. Whiting et al. 2003J, 25. ISchwartz et al 
20061 . 26. ISambruna et al.ll200ll , 27. 1 Jester et al.l 12002], 28. ISambruna et alj|2004 29. ISchwartz et alj|2000l . 30. IChartas et al 



2000 



Fig. IHb) plots the values of (Nx~d)max (= D 2 /d^ in , 
for convenience) as a function of the deprojected length of 
jets L. Interpolating them yields the scaling (A?A~d)max ~ 
10 9 L 9 p 6 c , indicating the trend that the capacity of outer- 
scale filaments increases as L increases (shaded area). In 
order to see the significance, we recall a promising scenario 
in which jets having large-scale magnetic fields are neces- 
sarily accompanied by hu ge currents in the bulk plasmas 
l|Appl fc Camenzindl Il992l ). the kinetic energy of which is 
dissipated o nly a little during the transport from the cores 
to the lobes l|Tashiro fc I sobc 2004). As a whole, this is com- 
patible with the superconductivity of plasmas. 

For example, let us consider the well-confirmed sam- 
ples M87 and Cen A, the nuclei of which have the super- 
mass ive black holes with mass ~ 10 9 M^ dMacchetto et all 
Il997l ) and (1O 7 -1O 8 )M dMarconi et al ll200rj). respectively . 
According to the arguments of lAppl fc Camenzindl (|l992T l. 
their central engines, incorporated with their accretion discs, 
ought to have the potential to drive a current of the order 
of magnitude of I ~ 10 19 A and ~ 10 18 A, respectively. 
Such a huge current could not be transported by a sin- 
;le uniform column, on account of the current inhibition 
Hondall2007l ) . One possible solution is to allow for the pres- 
ence of many filam ents that eac h carry a current limited by 
io ~ (mc 3 /e)/9,Tj (lHondall2000l). It is note d that the value 
of io is independent of A ( Honda et al.|[2000l ). The number of 
current filaments can then be estimated as A r (~ //io) ~ 10 15 
and ~ 10 14 , respectively. Note that these values are just in 



the expected ranges of (N\^d)ma,x for M87 and Cen A (cf. 
Fig. |4]d). The outcome suggests that, if the actual capacity 
of the outer-scale filaments reaches the level ~ (N\^d)max, 
the filament cluster, which consists of the smaller filaments 
with size A < d, would not be closely packed in the jet. 
This seems to be qualitatively amenable to the observed ap- 
pear arice_olM ; he_non^^ the jets of 
M87 (IBiretta et al lll99ll; IMarshall et al.ll2002l ) and Cen A 
l|Kraft et alj|2002l ; lHardcastle et al.ll2007h . 

As for the BLLac object Mrk421, t he black hole m ass 
of which is considered to be ~ 10 6 Mq dXie et al.lll998l ) or 
more, one anticipates a current of / > 10 17 A. Regard- 
ing the strongly beaming flow with a narrow viewing an- 
gle in which 6 ~ Tj, we find the order of io ~ 10 4 <5 A (in 
the regime, also note the scaling v c ~ b(io/e) in equation 
[HJ. For 5 ^ 100, we accordingly have N > 10 11 , which is 
much larger than (A ?- A~d)max ~ 10 8 (Fig. |4p). The relation 
N » (A r A~d)max, which is in contrast to N ~ (A r A~d)max for 
the aforementioned FRI sources, can also be found in other 
BLLacs (Mrk501, for example). It is thus ensured that the 
jets accompanying these compact objects have many fila- 
ments smaller than d (and thereby ; smalle r than D), consis- 
tent with the assumption given in iHondal |2008). 



5 DISCUSSION AND CONCLUSIONS 

Another signature of the 'quantization' of magnetized cur- 
rent filaments can be speculatively found in Fig. [2] By in- 



Scaling of the synchrotron cut-off 7 




(a) 



<D E 

s 

3 * 

O > 

— 

>- O) 
O o 
.£ — 

E « 

= c 

c a> 

E E 
Es= 



(b) 



-3-2-1012345 
Size of emission region logD(pc) 



1 7 
1 5 

< 
i 

;i 3 
1 1 

9 
7 
5 











/ inner-scale filaments 
with the skin depth 








I 


FR I 


1 


— 1 ♦ - 


BL Lac 


1 


■ 




• 









-4 -3-2-10 1 2 
Deprojected length of jet log L(kpc) 



Figure 4. (a) The lower limit of the transverse correlation length 
of filamentary turbulence d m i n versus the emission-region size D. 
The shaded area between the spatial limit (dotted line) and a 
power-law scaling (solid line) indicates the allowed domain for 
the correlation length. A skin depth is also indicated (e.g. for a 
plasma density of 10 — 3 cm -3 ; dot-dashed line), (b) The maxi- 
mum number of outer-scale filaments (Nxr^d)max versus the de- 
projected length of the jet L. The shaded area (below the scal- 
ing; solid line) suggests the capacity of the outer-scale filaments. 
The dot-dashed line indicates an upper limit: D 2 divided by the 
square of the skin depth (for 10~ 3 cm -3 ). BL Lac: circles, bars; 
FRI: filled square, fat bars; GRQ: filled diamond; CDQ: open 
diamonds, dotted bars. 



voking an approximate scaling B oc D" 1 ^ 2 , we find that the 
column energy density of the magnetic field e;,, estimated as 
~ B 2 D/8ir, is nearly constant: 



66 I 



const. 



10 1 



ergs cm 



(14) 



Obviously, this is at odds with the intuitive prediction of 
et oc D~ 2 reflecting three-dimensional (3D) adiabatic ex- 
pansion, implying that 2D expansion, arguably in the trans- 
verse plane, is preferentially inhibited due to the pressure of 
the l arge-scale ma gnetic field that would act radially inward 
(e.g. iHondal 120091 ). At the same time, it is also suggested 
that the inhomogeneous magnetic fields in the interior of 
emission regions robustly sustain such an averaged energy- 
density level against the self-contraction pressure. Consid- 
ering this aspect, the degenerate nature reflected in equa- 
tion (|14[) might be ascribed to the reasoning that any mag- 
netized filaments can by no means overlap spatially (or oc- 
cupy the common space) so as to carry a summed current ex- 
ceeding the critical, quantized value of io [for Tj > 0(1)], to 



prevent unlimited collapse (Honda 2000; Ho nda et alJ feoOO). 
We note that this exclusive property could not be deduced 
from the classical hydrodynamic description of current co- 
alescence, according to which an archetypical J x B pinc h 
leads to self-similar collapse (e.g. IZhdanov fc Vlasovlll998l ). 
The analogy to the updated notion can be found in the 
real mechanism for the Fermi pressure in de generate matter, 
based on the Pauli exclusion principle (e.g. ChandrasekharJ 
Il967h . With this insight, the B — D property revealed in 
Fig. [2] can be regarded as a macroscopic appearance of the 
quantization in relativistic flowing plasmas (Fig.[T]), and the 
beamed emission with p c (< io/e; Section Q as an acute 
appearance of collective quanta. The origin of the £(, value 
as such (in equation 1 14 p is another issue of note, but the 
details are beyond the scope of this paper. 

In conclusion, in the context of the filamentary jet 
model, we have analytically provided a new scaling of the 
cut-off frequency v c of non-diffuse synchrotron emissions 
(equation I10p . The arguments expanded here are based 
on the plausible assumption that small-amplitude magnetic 
fluctuations are superimposed on the local mean fields that 
permeate through current filaments. It has been demon- 
strated that the v c level could be much higher than the pre- 
viously suggested one (equation[T}, consistent with many X- 
ray measurements. The expected range of the present break 
frequency (equation|6} appears to be comparable to the pre- 
vious 'cut-off' (in the nomenclature of earlier literatures). In 
terms of the validity, we have provided extended discussions 
of the smearing effect on the synchrotron spectrum at v c , the 
inverse Comptonization effect and the ageing effect, which 
could involve the degenerative redshift of v . 

By combining the results obtained with key observa- 
tional data of AGN jets, we have found the property that 
the minimum correlation length of filamentary turbulence 
is larger for larger objects. In particular, we address the 
issue that the correlation data for compact BL Lacs and 
FRI radio sources are likely aligned in a common scaling, 
probably reflecting the spatiotemporal evolution of AGNs. 
This trend is amenable to the AGN unification scheme, in 
which the FR I sources are envisaged as misaligned BL Lacs. 
From translating the correlation scaling to the number 
scaling of filaments, we could infer how the filaments are 
packaged in the jet interior, and extract the consequence 
that the current carrying AGN jets must possess many 
filaments smaller than the entire emission-region size; 
e.g. for a pc-scale blazar jet the number of filaments 
is ~ 10 11 and more. The findings, including a constant 
column magnetic energy density, support the conjecture 
that AGN jets ubiquitously involve filamentary morphology. 



REFERENCES 

Appl S., Camenzind M., 1992, A&A, 256, 354 
Biermann P. L., Strittmatter P. A., 1987, ApJ, 322, 643 
Biretta J. A., Stern C. P., Harris D. E., 1991, AJ, 101, 1632 
Blasi P., 2002, Astroparticle Phys., 16, 429 
Blazejowski M. et al., 2005, ApJ, 630, 130 
Brunetti G., Mack K.-H., Prieto M. A., Varano S., 2003, 
MNRAS, 345, L40 
Celotti A., Ghisellini G., 2008, MNRAS, 385, 283 



8 M. Honda 



Chandrasekhar S., 1967, An Introduction to the Study of 

Stellar Structure. Dover, New York 
Chartas G. et al., 2000, ApJ, 542, 655 
Cheung C. C, Wardle J. F. C, Chen T., 2005, ApJ, 628, 

104 

Dominici T. P., Abraham Z., Teixeira R., Benevides- 

Soares P., 2004, ApJ, 128, 47 
Gaisser T. K., 1990, Cosmic Rays and Particle Physics. 

Cabbridge Univ. Press, New York 
Hardcastle M. J., Birkinshaw M., Worrall D. M., 2001, MN- 

RAS, 326, 1499 
Hardcastle M. J., Worrall D. M., Birkinshaw M., 

Laing R. A., Bridle A. H., 2002, MNRAS, 334, 182 
Hardcastle M. J., Kraft R. P., Worrall D. M., 2006, MN- 
RAS, 368, L15 
Hardcastle M. J. et al., 2007, ApJ, 670, L81 
Harris D. E., Krawczynski H., 2006, ARA&A, 44, 463 
Harris D. E., Finoguenov A., Bridle A. H., Hardcas- 
tle M. J., Laing R. A., 2002, ApJ, 580, 110 
Hillas A. M., 1984, ARA&A, 22, 425 
Honda M., 2000, Phys. Plasmas, 7, 1606 
Honda M., 2007, MNRAS, 376, 871 
Honda M., 2008, ApJ, 675, L61 
Honda M., 2009, ApJ, 706, 1517 
Honda M., Honda Y. S., 2004, ApJ, 617, L37 
Honda M., Honda Y. S., 2005a, ApJ, 633, 733 
Honda Y. S., Honda M., 2005b, MNRAS, 362, 833 
Honda M., Honda Y. S., 2007, ApJ, 654, 885 
Honda M., Meyer-ter-Vehn J., Pukhov A., 2000, 

Phys. Plasmas, 7, 1302 
Jester S., Roser H.-J., Meisenheimer K., Perley R., 2002, 

A&A, 385, L27 
Kang H., Jones T. W., 2007, Astroparticle Phys., 28, 232 
Kataoka J., Stawarz L., 2005, ApJ, 622, 797 
Kataoka J., Takahashi T., Makino F., Inoue S., Made- 
jski G. M., Tashiro M., Urry C. M., Kubo H., 2000, ApJ, 
528, 243 

Kataoka J., Stawarz L., Aharonian F., Takahara F., Os- 

trowski M., Edwards P. G., 2006, ApJ, 641, 158 
Keel W. C, 1988, ApJ, 329, 532 

Kneiske T. M., Bretz T., Mannheim K., Hartmann D. H., 
2004, A&A, 413, 807 

Konopelko A., Mastichiadis A., Kirk J., de Jager O. C, 
Stecker F. W., 2003, ApJ, 597, 851 

Kraft R. P., Forman W. R., Jones C, Murray S. S., Hard- 
castle M. J., Worrall D. M., 2002, ApJ, 569, 54 

Krawczynski H. et al., 2004, ApJ, 610, 151 

Lobanov A. P., Zensus J. A., 2001, Sci, 294, 128 

Longair M. S., 1994, High Energy Astrophysics, Vol. 2: 
Stars, the Galaxy and the Interstellar Medium. Cam- 
bridge Univ. Press, Cambridge 

Macchetto F., Marconi A., Axon D. J., Capetti A., 
Sparks W., Crane P., 1997, ApJ, 489, 579 

Mack K.-H., Prieto M. A., Brunetti G., Orienti M., 2009, 
MNRAS, 392, 705 

Marconi A., Pastorini G., Pacini F., Axon D. J., Capetti A., 
Macchetto D., Koekemoer A. M., Schreier E. J., 2006, 
A&A, 448, 921 

Marshall H. L., Miller B. P., Davis D. S., Perlman E. S., 
Wise M., Canizares C. R., Harris D. E., 2002, ApJ, 564, 
683 

Marshall H. L. et al., 2010, ApJ, 714, L213 



Meisenheimer K., Roser H.-J., Schlotelburg M., 1996, 
A&A, 307, 61 

Meisenheimer K., Yates M. G., Roser H.-J., 1997, A&A, 
325, 57 

Montgomery D., Liu C. S., 1979, Phys. Fluids, 22, 866 
Miicke A., Protheroe R. J., 2001, Astroparticle Phys., 15, 
121 

Ostrowski M., 1988, MNRAS, 233, 257 
Owen F. N., Hardee P. E., Cornwell T. J., 1989, ApJ, 340, 
698 

Perez-Fournon I., Colina L., Gonzalez-Serrano J. I., Bier- 
mann P. L., 1988, ApJ, 329, L81 

Perlman E. S., Biretta J. A., Sparks W. B., Mac- 
chetto F. D., Leahy J. P., 2001, ApJ, 551, 206 

Pesce J. E., Sambruna R. M., Tavecchio F., Maraschi L., 
Cheung C. C, Urry C. M., Scarpa R., 2001, ApJ, 556, 
L79 

Sambruna R. M., Urry C. M., Tavecchio F., Maraschi L., 
Scarpa R., Chartas G., Muxlow T., 2001, ApJ, 549, L161 

Sambruna R. M., Gambill J. K., Maraschi L. M., Tavec- 
chio F., Cerutti R., Cheung C. C, Urry C. M., Char- 
tas G. C, 2004, ApJ, 608, 698 

Schlickeiser R., 2002, Cosmic Ray Astrophysics. Springer, 
Berlin 

Schwartz D. A. et al., 2000, ApJ, 540, 69 

Schwartz D. A. et al., 2006, ApJ, 640, 592 

Stawarz L., Sikora M, Ostrowski M., 2003, ApJ, 597, 186 

Tansley D., Birkinshaw M., Hardcastle M. J., Wor- 
rall D. M., 2000, MNRAS, 317, 623 

Tashiro M., Isobe N., 2004, Astron. Herald, 97, 400 

Tavecchio F. et al., 2001, ApJ, 554, 725 

Tingay S. J., Lenc E., Brunetti G., Bondi M., 2008, ApJ, 
136, 2473 

Urry C. M., Padovani P., 1995, PASP, 107, 803 
Urry C. M. et al., 1997, ApJ, 486, 799 
Waters C. Z., Zepf S. E., 2005, ApJ, 624, 656 
Whiting M. T., Majewski P., Webster R. L., 2003, PASA, 
20, 196 

Worrall D. M., Birkinshaw M., 2005, MNRAS, 360, 926 
Xie G. Z., Bai J. M., Zhang X., Fan J. H., 1998, A&A, 334, 
L29 

Zhang Y. H. et al., 2002, ApJ, 572, 762 
Zhdanov S. K., Vlasov V. P., 1998, Plasma Phys. Rep., 
24, 497 



This paper has been typeset from a TpX/ DTjjX file prepared 
by the author. 



